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0.1 Abstract 

The Jaynes-Cummings model of a two-level atom in a single mode cavity is of 
fundamental importance both in quantum optics and in quantum physics gen- 
erally, involving the interaction of two simple quantum systems - one fcrmionic 
system (the TLA), the other bosonic (the cavity mode). Depending on the 
initial conditions a variety of interesting effects occur, ranging from ongoing 
oscillations of the atomic population difference at the Rabi frequency when the 
atom is excited and the cavity is in an n-photon Fock state, to collapses and 
revivals of these oscillations starting with the atom unexcited and the cavity 
mode in a coherent state. The observation of revivals for Rydberg atoms in a 
high-Q microwave cavity is key experimental evidence for quantization of the 
EM field. Theoretical treatments of the Jaynes-Cummings model based on ex- 
panding the state vector in terms on products of atomic and n-photon states 
and deriving coupled equations for the amplitudes are a well-known and simple 
method for determining the effects. 

In quantum optics however, the behaviour of the bosonic quantum EM field 
is often treated using phase space methods, where the bosonic mode annihila- 
tion and creation operators are represented by c-number phase space variables, 
with the density operator represented by a distribution function of these vari- 
ables. Fokker-Planck equations for the distribution function are obtained, and 
either used directly to determine quantities of experimental interest or used 
to develop c-number Langevin equations for stochastic versions of the phase 
space variables from which experimental quantities are obtained as stochastic 
averages. Phase space methods have also been developed to include atomic 
systems, with the atomic spin operators being represented by c-number phase 
space variables, and distribution functions involving these variables and those 
for any bosonic modes being shown to satisfy Fokker-Planck equations from 
which c-number Langevin equations are often developed. However, atomic spin 
operators satisfy the standard angular momentum commutation rules rather 
than the commutation rules for bosonic annihilation and creation operators, 
and are in fact second order combinations of fermion annihilation and creation 
operators. Though phase space methods in which the fermion operators are 
represented directly by c-number phase space variables have not been success- 
ful, the anti-commutation rules for these operators suggests the possibility of 
using Grassmann variables - which have similar anti-commutation properties. 
However, in spite of the seminal work by Cahill and Glauber and a few appli- 
cations, the use of phase space methods in quantum optics to treat fermionic 
systems by representing fermionic annihilation and creation operators directly 
by Grassmann phase space variables is rather rare. This paper shows that phase 
space methods using a positive P type distribution function involving both c- 
number variables (for the cavity mode) and Grassmann variables (for the TLA) 
can be used to treat the Jaynes-Cummings model. Although it is a Grassmann 
function, the distribution function is equivalent to six c-number functions of 
the two bosonic variables. Experimental quantities are given as bosonic phase 
space integrals involving the six functions. A Fokker-Planck equation involving 



both left and right Grassmann differentiation can be obtained for the distribu- 
tion function, and is equivalent to six coupled equations for the six c-number 
functions. 

The approach used involves choosing the canonical form of the (non-unique) 
positive P distribution function, in which the correspondence rules for the 
bosonic operators are non-standard and hence the Fokker-Planck equation is 
also unusual. Initial conditions, such as those above for initially uncorrelated 
states, are discussed and used to determine the initial distribution function. 
Transformations to new bosonic variables rotating at the cavity frequency en- 
ables the six coupled equations for the new c-number functions - that are also 
equivalent to the canonical Grassmann distribution function - to be solved an- 
alytically, based on an ansatz from an earlier paper by Stenholm. It is then 
shown that the distribution function is exactly the same as that determined 
from the well-known solution based on coupled amplitude equations. 

In quantum-atom optics theories for many atom bosonic and fermionic sys- 
tems are needed. With large atom numbers, treatments must often take into 
account many quantum modes - especially for fermions. Generalisations of 
phase space distribution functions of phase space variables for a few modes to 
phase space distribution functionals of field functions (which represent the field 
operators, c-number fields for bosons, Grassmann fields for fermions) are now 
being developed for large systems. For the fermionic case, the treatment of the 
simple two mode problem represented by the Jaynes-Cummings model is a use- 
ful test case for the future development of phase space Grassmann distribution 
functional methods for fermionic applications in quantum-atom optics. 



1 Introduction 

The simplest theoretical model in quantum optics involves a two state atomic 
system coupled to a single mode quantum electromagmetic (EM) field, in which 
the system transition frequency is in near resonance with the mode frequency. 
This system was first treated using the rotating wave approximation and is re- 
ferred to as the Jaynes-Cummings model p. The model can be used not only 
in cavity quantum electrodynamics (QED) to describe a two level atom (TLA) 
in a high-Q cavity - which could be two Rydberg atom states in a microwave 
cavity, or two lower lying atomic states in an optical cavity, but also in laser 
spectroscopy for treating a two level atom in free space coupled to a single mode 
laser field in a coherent state. The Jaynes-Cummings model is of fundamental 
importance both in quantum optics and in quantum physics generally, involv- 
ing the interaction of two simple quantum systems - one being fermionic (the 
TLA) , the other bosonic (the cavity mode) . In laser spectroscopy the prediction 
[2], interpretation [3] and observation [4], [5], [6] of the three peak resonance 
fluorescence spectrum for a TLA strongly driven by a narrow bandwidth laser 
field stimulated much interest in the new research field of quantum optics. In 
the cavity QED case, a variety of interesting effects occur depending on the 
initial conditions, ranging from ongoing oscillations of the atomic population 
difference at the Rabi frequency when the atom is excited and the cavity is in 
an n-photon Fock state, to collapses and revivals of these oscillations starting 
with the atom unexcited and the cavity mode in a coherent state [7], [5], [S]. 
The collapse time corresponds to the time scale for phase factors with photon 
numbers at the extremes of the photon number distribution to get out of phase 
by TT. The revival time corresponds to when the phase factors for neighboring 
photon numbers get out of phase by 27r, and hence rephase. For the coherent 
state case the revival time is \f^ times longer than the collapse time, where n 
is the mean photon number of the coherent field. The observation of revivals 
[lOj . |llj for Rydberg atoms in a high-Q microwave cavity is key experimental 
evidence for quantization of the EM field. Over 1200 papers have been writ- 
ten about the Jaynes-Cummings model indicating its importance in quantum 
optics, and a comprehensive review is provided by Shore and Knight |12j . 

Theoretical treatments of the Jaynes-Cummings model based on expanding 
the state vector in terms of products of atomic and n-photon states and deriving 
coupled equations for the coefficients are a well-known and simple method for 
determining the effects [13] , [9] . Another simple approach is based on using the 
actual atom-cavity mode energy eigenstates - the dressed states [3] - to describe 
the behaviour. In Stenholm's approach |14j the state vector was expanded in 
products of Bargmann states for the field and spinor states for the TLA. How- 
ever, in quantum optics the behaviour of the bosonic quantum EM field is often 
treated using phase space methods [15], [16], [13], [17] with the density oper- 
ator represented by a phase space distribution function and the bosonic mode 



annihilation and creation operators represented by c-number phase space vari- 
ables. Fokker-Planck equations for the distribution function are obtained, and 
either used directly to determine quantities of experimental interest or used to 
develop c-number Langevin equations for stochastic versions of the phase space 
variables, such as when distribution functions of the positive P type [IH] are 
used. Phase space methods have also been developed to include atomic systems 
[15] ■ |16) . with the atomic spin operators being represented by c-number phase 
space variables, and distribution functions involving these variables and those for 
any bosonic modes being shown to satisfy Fokker-Planck equations from which 
c-number Langevin equations are often developed. However, atomic spin opera- 
tors satisfy the standard angular momentum commutation rules rather than the 
commutation rules for bosonic annihilation and creation operators, and are in 
fact second order combinations of fermion annihilation and creation operators. 
It would perhaps be more natural to use a phase space method in which the 
phase space variables represent the fermion operators directly. The problem 
however is that it is not possible for fermion annihilation and creation operators 
to be represented by c-number variables. The anti-commutation rules for these 
operators suggests the use of variables which have similar anti-commutation 
properties, and such variables are Grassmann numbers. These numbers were 
invented in the 19th century by the mathematician Hermann Grassmann, but 
their properties would be unfamiliar to many physicists. A general account of 
Grassmann variables is to be found in the book by Berezin [19] . 

Apart from the seminal work by Cahill and Glauber |20) and a few appli- 
cations [21], [22], [23], [24] the use of phase space methods in quantum optics 
to treat fermionic systems by representing fermionic annihilation and creation 
operators with Grassmann variables [1^, [20] is rather rare. In contrast, this 
approach is widely used in particle physics [IS] , [33] , [37] . This paper shows how 
phase space methods using a distribution function involving both c-number vari- 
ables (for the cavity mode) and Grassmann variables (for the TLA) can be used 
to treat the Jaynes-Cummings model. Quantities of experimental interest such 
as atomic state populations and coherences or mean numbers of photons are 
related to quantum correlation Junctions, which are expectation values of nor- 
mally ordered products of boson, fermion annihhilation and creation operators. 
Accordingly, the distribution function chosen is of the positive P type for the 
bosonic and is similar to the complex P type for the fermionic variables (for 
short, a positive P distribution). Unlike the more usual Glauber-Sudarshan P, 
Husimi Q or Wigner W distributions, the positive and complex P distributions 
involve a double phase space, with two phase space variables for each bosonic or 
fermonic mode. The distribution function is not unique but the existence of the 
distribution function in a canonical form can be established. The distribution 
function is a Grassmann function, and is specified by six c-number functions of 
the two bosonic variables. Experimental quantities are given as bosonic phase 
space integrals involving the six functions. The non-uniqueness of the positive 
P distribution is reflected in different Fokker-Planck equations based on differ- 
ing correspondence rules. Fokker-Planck equations involving both left and right 
Grassmann differentiation can be obtained for both the canonical and stan- 



dard distribution functions, and these are equivalent to six coupled equations 
for the six c-number bosonic functions. Initial conditions, such as for initially 
uncorrelated states, are discussed and used to determine the initial canonical 
distribution function. The positive P type distribution functions would en- 
able the phase space variables to be replaced by stochastic variables and the 
Fokker-Planck equation replaced by Langevin stochastic equations of the Ito 
type. However, this development is not needed for the present paper as the 
Fokker-Planck equation for the canonical distribution function is solved directly 
via an adaption of Stenholm's method [13] . The solution is entirely in agreement 
with the standard quantum optics results. However, the analogous solution of 
the Fokker-Planck equation for the standard distribution function leads to a dis- 
tribution function that diverges on the phase space boundary and for large times, 
indicating that the original Fokker-Planck equation was not valid. Applications 
of the positive P distribution using the canonical distribution function may be 
preferable to those where the derivation of the Fokker-Planck equation is based 
on the standard correspondence rules. For the canonical distribution function 
more general Fokker-Planck equations involving derivatives higher than second 
order may occur [28], so that no replacement by Langevin stochastic equations 
is possible. The Fokker-Planck equation may also not have a positive definite 
diffusion matrix. 

The use of phase space approaches to the Jaynes-Cummings model is itself 
rather unusual, however the paper by Eiselt and Risken ^S] does use such a 
method. A generalised Wigner distribution for the field mode is constructed 
from four field mode operators given by the matrix elements of the full density 
operator with the two atomic states, then combined with the field mode an- 
nihilation and creation operators. A generalised density operator for the field 
with six components is obtained and used to define a six component character- 
istic function and hence a six component Wigner function. Damping effects are 
also included. For zero damping, collapse and revival effects are found for the 
population difference of the atomic states. We note that their approach also 
requires consideration of six field distribution functions, as is required in the 
present work. However, their approach cannot be related in any simple way to 
ours. 

Although it is of interest to show that the Jaynes-Cummings model can be 
treated via Grassmann phase space methods the work presented here has more 
general relevance. In quantum-atom optics Bose-Einstein condensates and de- 
generate Fermi gases are being studied and theories for many atom bosonic and 
fermionic systems are needed. With large atom numbers, treatments must of- 
ten take into account many quantum modes - especially for fermions in view of 
the Pauli exclusion principle. Phase space distribution function methods based 
on separate variables for each mode then start to become unwieldy. Gener- 
alisations of phase space distribution function methods involving phase space 
variables for a few modes to phase space distribution functional methods based 
on field functions which represent the entire field operators are now being devel- 
oped for large systems. The field functions are c-number fields for the bosonic 
case, Grassmann fields for the case of fermions. This development is more ad- 



vanced for bosonic systems, with applications already made to problems such 
as optical light in dispersive non-linear media [30], [31], spatial squeezing in 
quantum EM fields [35] and quantum noise in Bose-Einstein condensates [55] . 
However, distribution functional methods involving atomic media have also been 
formulated |34) . based on a different approach in which the atomic system is 
treated in terms of c-number fields representing the atomic spin field opera- 
tors rather than the particle field creation and annihilation operators. More 
recently, an alternative phase space approach for fermions, the Gaussian oper- 
ator basis method [3S], originally developed for bosonic systems [33] has been 
published. Here the fermionic density operators are represented as a positive 
distribution over a generalised phase space. The phase space variables again are 
c- numbers, now associated with pairs of fermion annihilation and/or creation 
operators, and Grassmann numbers are only used to establish properties of the 
Gaussian operators. However, less progress has been made for the different ap- 
proach - analogous to that applied to bosonic systems - in which the fermion 
particle field operators are represented by Grassmann fields, and the treatment 
of the simple two mode problem represented by the Jaynes-Cummings model 
is a useful first step in the development of this type of phase space Grass- 
mann distribution functional methods for applications in quantum-atom optics. 
For such fermion multimode systems, functional Fokker-Planck equations would 
be obtained for the distribution functional, and these would be converted into 
equivalent Langevin stochastic field equations of the Ito type. It is expected 
that the latter step would be important for practical numerical applications, 
since solving for the full distribution functional would likely be impractical, and 
unnecessary if the required physical predictions can be made via less computer 
intensive stochastic methods. 

In section [5] the basic physics of the Jaynes-Cummings model is reviewed, 
with the Hamiltonian, the physical and spin states, atomic population and tran- 
sition operators, photon number operators and spin operators being treated in 
terms of fermion and boson creation, annihilation operators. The phase space 
approach is developed in section [S] starting with defining the quantum correla- 
tion functions then introducing the characteristic function and its phase space 
integral relation to the distribution function. The characteristic and distribution 
functions are shown to be certain Grassmann functions with c-number coeffi- 
cients that are functions of the bosonic phase space variables, and the expression 
for the canonical form of the distribution function is presented. Results for the 
quantum correlation functions, atomic populations and coherences are obtained 
in terms of phase space integrals. In section|3]correspondence rules for the effects 
of boson and fermion operators on the density operator are obtained, and used 
to derive Fokker-Planck equations for the Jaynes-Cummings model. The stan- 
dard and canonical correspondence rules are presented and the Fokker-Planck 
equation is obtained for the canonical distribution function. Sets of coupled 
equations for the c-number distribution function coefhcients are obtained from 
the Fokker-Planck equations. The original c-number distribution function co- 
efficients are replaced by coefficients in a frame rotating at the cavity freqency 
and expressed in terms of rotating phase space variables. The initial form of 



the distribution function is found for the case of uncorrelated initial states, such 
as those leading to effects such as Rabi oscillations and collapse/revival phe- 
nomena. In section [5] the coupled c-number equations are solved analytically 
for the canonical distribution case and shown to predict the characteristic the 
Jaynes-Cummings mode effects. The conclusions for the paper are set out in 
section [51 Basic results for Grassmann algebra and calculus and for boson and 
fermion Bargmann states are presented in the Appendices. 

2 Basic Physics of Jaynes-Cuminings Model 

2.1 Atomic System and Cavity Mode 

The basic atomic system treated in the Jaynes-Cummings model is a single 
atom system with just two internal states, denoted |1), |2), with energies -Ei, 
E2. Centre of mass degrees of freedom are ignored. There are therefore just four 
atomic operators to consider, denoted ct,;j = \i) {j\. Two of these are population 
operators. Pi = |1) (1| and P2 ~ |2) (2|, and two are atomic transition operators 
(7+ = |2) (1| and (t_ = |1) (2|. In the simple Jaynes-Cummings model the two- 
level atom interacts with a cavity mode. Let a, a' be the boson annihilation, 
creation operators for the field mode. The states for the cavity mode are the 
usual n photon states, denoted \n), with n = 0, 1, 2, ... If the atomic transition 
frequency is ujq = (£'2 — Ei)/h, the frequency for the field mode is w and the 
atom-field mode coupling constant (one-photon Rabi frequency) is fl, then the 
Hamiltonian for the one atom Jaynes-Cummings model is given by 

Hjc = EAiP2 +Pi) + ^hujoiP2 - Pi) + huj{a^a) + hn{d+a + a^d_) (1) 

where the average atomic energy is Ea — {E2 + Ei)/2. The first term is usually 
ignored as it only introduces a phase factor exp(— i?^t/?i) into the evolution of 
all states in the one atom Jaynes-Cummimgs model. 

The dynamics of the one atom Jaynes-Cummings model is treated in many 
quantum optics textbooks and papers (see ^7\, [T3], [5]). The state vector is 
expanded in terms of n photon states for the field and the two internal atomic 
states as 

L*) = exp{-EAt/h) (2) 

X Y^ i ^;rrrT2(^) cxp(-i(n - Iw + -ojo)t) |2) |n - 1) + A„i{t) exp{~i{nuj - -u}o)t) |1) \n) j 

involving interaction picture amplitudes A„i , ^;jirf2 • '^^^'^ coupled equations for 
the amplitudes 

d 1 

«^^;r3T2 = -il^/ncxp{+iAt) A^i 

d 1 

j^^ni = -^J^/r^cxp(-^Ai) A:;—^^ (3) 



involve the detuning 

A = ujQ- u} (4) 

and may be solved using Laplace transforms. The solution is 

1/1 1 r A A„i(0)-f2V7IA;^3Y2(0) 

Anl{t) = cxp{--iAt) [coa{-uJnt){Ani{0)} +lSm{-UJnt) 



n^ Ani{0) + A A^^^iO) 



-n-i2yj - exp(+-iAf) fcos(-a;„t){A;j^2(0)}-«sin(-a;„t) ^ ^ 



where 



v/A2 + nil^ (6) 



is the Rabi frequency. Population and coherence oscillations at various frequen- 
cies ijJn are therefore predicted, and which frequencies are observed depends on 
the inital conditions. For example, if the atom is not excited and the field is 
in an m photon state, then A„i(0) = Sum, ^71312(0) = ^-^d hence oscillations 
at a single frequency uJm will occur. On the other hand, if the field was in a 
coherent state with amplitude 77, then oscillation frequencies clustered around 
the mean frequency uJn associated with mean photon number n = \'i]\'^ occur, 
the range of frequencies being associated with the standard deviation An = v^ 
in photon numbers. This results in the collapse, revival effects discussed previ- 
ously. In an alternative approach, Stenholm '14J expanded the state vector in 
products of Bargman states for the field and spinor states for the TLA and also 
demonstrated collapse, revival effects. As we see, his solution can be adapted 
to determining the Grassmann distribution function in the present treatment. 



2.2 Fermion and Bosonic Modes 

Instead of treating the Jaynes-Cummings model via the usual elementary quan- 
tum optics approach - for example via matrix mechanics with basis states |1) |n), 
|2) \n — 1), we can replace the original one two-level atom plus cavity mode sys- 
tem with a somewhat enlarged system consisting of two fermionic modes 1, 2 
interacting with one bosonic mode. This enlarged system includes states that 
are in one-one correspondence with the states for the original one two-level 
atom plus single cavity mode system, and the general dynamics for the larger 
Fermi-Bose system incorporates all possible behaviour for the original Jaynes- 
Cummings model. Note that relating the two state atomic system to two fermion 
modes has nothing to do with whether the atom itself is a fermion or a boson. 
This process involves the introduction of fermion annihilation, creation op- 
erators Ci,cl {i — 1,2) for the two fermion modes, which satisfy standard fermion 
anti-commutation rules. The commutation and anti-commutation rules for the 



(5) 



ferniion and boson operators are 



{ct,cl} = S 



3 J - "«j 



{q,c,} = {cT,c}-0 



-t ^ti 



[a,a1 = 1 (7) 

with the boson and fermion operators commuting. 

We then introduce the four possible fermion Fock states 

|mi;m2) = (cIn(ctrM0> (8) 

where mi, m2 are equal to 0, 1 only and |0) is the vacuum state. For the vacuum 
state Ci |0) =0, (0| c| ~ 0. Thus |0; 0) will be the vacuum state, with no fermions 
in either mode, |1; 0) and |0; 1) will be the one fermion states with one fermion 
in modes 1,2 respectively, and |1;1) will be the two fermion state with one 
fermion in each mode. Naturally, the Pauli exclusion principle excludes states 
with more than one fermion in any mode. 

The full basis states for the Jaynes-Cummings model are Fock states of the 
form 

|mi;m2;n) = (clrH4r^^|0) (9) 



where mi = 0, 1 give the number of atoms in state 2 = 1,2 and n = 0, 1, 2, ... 
gives the number of photons in the field mode. Obviously since (cj)^ = there 
can not be more than one atom in state i. The state with no atoms or photons 
present is the vacuum state |0). 



2.3 Quantum States 

The most general mixed state for which there may be either zero, one or two 
atoms has a density operator 

P = 5I(^oo«;OOm|0;0;n) (0;0,m| 

n.m 

+ J2^Pwn;iOm U; 0; n) (1; 0, m\ + Pio„;oim |1; 0; "> (0; 1, "i| 

n,m 
+P01n;10m \^i 1; n) (1; 0, to| + Poin-film l^', 1; «) (0; 1, m|) 
+ J2(Piin;llm\^-^;n){l;l,m\ (10) 

where the Pijn-kim (hj^kj = 0,1) are density matrix elements. This corre- 
sponds to a statistical mixture of states in which there are no atoms, one atom 
or two atoms. Pure states which are quantum superpositions of states with 
differering numbers of atoms are forbidden by the super-selection rule based 
on conservation of atom number. Statistical mixtures of such forbidden states 
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would in general lead to the presence of coherences between states with differing 
numbers of atoms. Thus there are no coherences between states with different 
numbers of atoms in (ITU|) . 

The condition that the trace of the density operator is one requires 

Trp — 2_^(P00n;00n + PlOn;10ra + Poira;01n + Plln;lln) = 1 (H) 

n 

Note that in addition to atomic states corresponding to the one atom Jaynes- 
Cummings model, we now have two extra states in our enlarged system, and 
we need to consider whether these have any physical significance. The vacuum 
state |0; 0) could correspond to an experiment with only the cavity field present 
but with no atoms inside the cavity - not a situation of great interest. The two 
fermion state |1; 1) could correspond to a two atom system inside the cavity, 
with one atom in state |1) and the other in state |2). All four fermion states 
\'mi]in2) also have a mathematical interpretation in terms of spin states in our 
enlarged system. 



2.4 Normal Ordering and Vacuum State Projector 

The vacuum state |0) which contains no bosons or fermions in any mode is 
associated with the vacuum state projector |0) (0|. It turns out that the vacuum 
projector can be written in terms of normally ordered forms of products of 
exponential operators based on mode number operators. The normally ordered 
form for a product of annihilation, creation operators involves placing all creation 
operators to the left of all annihilation operators, whilst retaining all creation 
and all annihilation operators in their original order, and multiplying by either 
+1 or —1 depending on whether the required permutation oi fermion operators 
from the original product is even or odd .38;. 

General functions of the annihilation, creation operators (such as the Hamil- 
tonian) can always be written as linear combinations of products with c-number 
coefficients, and the normally ordered form of such a function is just the linear 
combination of the normally ordered form of the various products. A useful the- 
orem that applies to the normally ordered product of two functions involving 
the number operators for two different fermion modes is 

N{F{c^fi,) G{Z]Z,)) = Af{F{clc,))Af{G{d]d,)) [i ^ j) (12) 

where F, G are arbitrary functions. This result only depends on the annihilation 
and creation operators anticommuting for i ^ j. 

The vacuum state projector is given in terms of normally ordered operators 
as 

|0)(0| = AA(nexp(-aIa,))=AA(cxp(-^a|a,)) (13) 

i i 

= AA(nexp(-clQ))=AA(exp(-^c|B,)) (14) 
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for bosons and fermions respectively. The ([T^ result for bosons is given in [T3] 
(see p. 45). The result for fermions follows from the vacuum state requirements 
that Ci |0) — 0, (0|cj = are also satisfied by the expression in (ITU) , together 
with the requirement that the vacuum state matrix element for the projector is 
also unity. 

2.5 Population and Transition Operators 

We now identify the two one fermion Fock states with the two internal atomic 
states in the one atom Jaynes-Cummings model as follows: 

cl|0> = |1;0>^|1) 

cUo) = |0;1>^|2) (15) 

and identify the four different products involving an annihilation, a creation op- 
erator and the vacuum state projector with the four one atom atomic population 
or projector operator Pi 2 and transtion operators a± 

d\ |o) (o| di ^=^ Pi cl |o) (o| C2 ^=^ P2 

d\ |0) (0| ?2 ^=^ d- 4 |0) (0| ?i ^=^ d+ (16) 

For the expanded Jaynes-Cummings model the two atom state population or 
projector operator P12 with one atom in each of the atomic states we have 

c\cl\0) = \l;l) (17) 

ctcj |0) (0| C2?i ^^ P12 (18) 

and for the zero atom population or projector operator Pq with no atoms in 
either of the internal atomic states we have 

|0) = |0;0) (19) 

|0) (0| ^=^ Po (20) 

This of course is the same as the vacuum state projector. 

Using Eq. ^^ and the result P^ we can obtain expressions for all the atom 
state operators just in terms of fermion annihilation, creation operators. We 
find for the one atom population and transition operators 

Pi = CJCI - c{c2C2Ci P2 = clc2 ~ c\c\c2Ci (21) 

(T_ = c\c2 (?+ = c|ci (22) 

and for the two atom and zero atom population operators 

Pi2 - c\Zlc2Ci (23) 

Po = 1 - c\ci — c\c2 + c\c\c2Ci (24) 
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The population operators thus satisfy the condition 

P0 + P1+P2 + P12 = 1 (25) 

which (as wih be seen below), just represents the total probability for any general 
mixed state (jlOl) that the system is found in either a zero, one or two atom state 
being equal to unity 

2.6 Hamiltonian and Number Operators 

The Hamiltonian needs to be modified for the enlarged system to allow for states 
with two or zero atoms. The one atom Jaynes-Cummings model Hamiltonian 
Hjc is replaced by an enlarged Hamiltonian H by adding a two atom energy 
term {E2 + Ei)Pi2 to the Hamiltonian in ([T|). This enlarged Hamiltonian has 
zero matrix elements between two atom states |1; 1), one atom states |0; 1), |1; 0) 
and zero atom states |0; 0). The two and zero atom states are eigenstates of the 
Hamiltonian with atomic energies E2 + Ei and 0, as expected. The atom-field 
coupling term also cannot cause a transtion between states with differing total 
atom number, consistent with not violating the super-selection rule. 

In terms of the bosonic and fermionic annihilation, creation operators the 
enlarged Hamiltonian for the Jaynes-Cummings model may be written as 

H — Ejl,{c2C2 + c\ci) + —hLjQ{c2C2 — c\c^ (26) 



where the results ((2T|) . ((22|) and ((23|) for the one and two atom operators have 
been substituted into the enlarged Hamiltonian. Terms involving c}c2C2Ci cancel 
out. 

The number of photons present is determined from the operator 



- n'in (27) 



n — a'a 



whose eigenvalues are n = 0, 1, 2, .We can also introduce number operators for 
the two atomic states via 

n^ = cld, (i = l,2) (28) 

where from the femion anticommutation rules the eigenvalues for the atomic 
number operators are 0, 1 only. 

The total number of atoms present is given by the number operator 

N = {clc2+c\ci) (29) 

= X Po + 1 X (Pi + P2) + 2 X Pi2 (30) 

the second expression for the number operator being related to the projectors 
Pq, Pi, P2 and P12 being for zero, one, one and two atom states respectively. 
The number operator can have eigenvalues 0,1,2 and clearly 

N \mi; m2; n) = {mi + m2) |mi;m2;n) (31) 
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The number operator commutes with the Hamiltonian, so if we initially have 
a physical state with a specified number of atoms, then state evolution does not 
change the atom number. Thus one atom states do not change into two atom 
states and if the initial density operator corresponds to the form in (ITUl) with 
PoonOOm = Piiniim = ^'^^ One atom states it will remain in this form. There 
is of course no conservation law for the photon number. . 

2.7 Probabilities and Coherences 

We can obtain expressions for the physical quantities in terms of the annihila- 
tion, creation operators. For the one atom probability of finding one atom in 
state |1) and none in state |2) we have 



Pi = Tr{Pi p) = (d\di) - (?I?^?2Ci) = Tr{dip8{) - Tr(?2?iP?l?^) (32) 



where ( ^) = Tr(p^). Similarly the one atom probability of finding one atom 
in state |2) and none in state |1) is 

P2 = Tr{P2 p) = (clc2) - (cIct?2Ci\ = Tr(?2P?J) - Tr(?2Cipclc|) (33) 






Note that these probabilities are not the same as (clcij as might be expected. 
The one atom coherences between state |1) and state |2) are given by 

P12 = Tr(B^p) = (c|c2) = Tr{d2pc\) 



P21 = Tr(d+ p) = (dldi) ^Tr{di^cl) (34) 






The mean number of photons is 

n = (a'' a) (35) 

Also, the two atom probability P12 for finding one atom in state |1) and one 
in state |2), and the zero atom probability Pq for finding no atom cither in state 
|1) or in state |2) are given by 



P12 = Tr(Pi2p) = (^c|c^C2Ci) = Tr(c2Cipclc5) (36) 

Po = TriPop) = l-Pi-P2-Pi2 (37) 

the last result just expressing the fact that all the probabilities must add up to 
one. Expressions involving the density matrix elements for general mixed states 
can be easily obtained for all these results. 

Clearly, for one atom states where PoonOOm = Piin-iim = we have using 

m 

Po = Pl2 = Pl+P2 = l (38) 

" = X!"(^10";10«+P01«;01n) (39) 
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For general mixed states as in (ITUl) we see that expectation values of odd 
numbers of fermionic creation and destruction operators are zero. Thus 

cl) = (cI4?,) = (cl?2Ci)=0 (40) 

For one atom states the only non-zero expectation values are those that involve 
one fermion annihilation operator and one creation operator of the form ( c-c-^ 



Thus in addition to the last results the expectation value involving four fermion 
operators is zero 

(clc^C2?l)-Pl2-0 (41) 

since Pun-iin ^ for one atom states. This corresponds to the probability of 
finding two atoms present being zero in the one atom Jaynes-Cummings model. 
For general mixed states the quantity (C1C2C2C1) is conserved, since CJC2C2C1 
commutes with the Hamiltonian 

Although we will be mainly focused on states corresponding to one atom 
states in the Jaynes-Cummings model, it will be convenient to introduce non- 
physical states involving coherent superpositions of different total numbers of 
atoms and even states where the expansion coefficients are Grassmann numbers. 
For non-physical states involving superpositions of different atom numbers there 
may be off diagonal elements of the density matrix between states of different 
atom numbers, and the results in Ea.(|40)) may not apply. 



3 Phase Space Theory 

There are several ways to set out phase space theory. Here we start with the 
characteristic function - which is equivalent to the set of all normally ordered 
quantum correlation functions - and then relate this to the distribution function 
- whose phase space integral with products of phase space variables determines 
the correlation functions. The normally ordered quantum correlation functions 
are defined as 

G(mi,m2,ri; p,l2,h) 

where mi,li = 0, 1 only. Normally ordered correlation functions in which the 
density operator appears in the middle of the trace with annihilation operators 
on the left and creation operators on the right appear in fundamental treatments 
of boson and fermion detection processes (see for example [39]). 
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3.1 Characteristic Function 

We define the characteristic function x(Cj? ; ^i; f^t) ^^^ 

x{^,e,h,h+) = Tiin+{e)^+{h+)pnj{h)n-{0) (43) 

n+{^+) = expzaC+ f2j7(C) =expi^a^ (44) 

n+{h+) = expi ^ Cj/i+ 17j(/i) = expi ^ /iic| (45) 

4=1,2 1=1,2 

For the bosonic mode we associate a pair of c-numhers ^, ^ . For each fermionic 
mode i = 1,2 we associate a pair of Grassmann-numhers hi,hj and h = 
{/ii, /12}, /i^ = {/ij^, /i^}. The characteristic function will be a c- number analytic 
function of ^,^ and a Grassmann junction of /ii, /i]^, /12, ^2 • The factors for 
the different modes in the last two expressions for r2't(/i+) and (ij {h) may be 
put in any order since they commute, but by convention the order will be 2, 1 
for rit(/i+) and ,1,2 for Clj{h). This leads to a quasi distribution function of 
the positive P type, which is generally better suited for obtaining Ito stochastic 
equations. 

For the physical states where the state vector involves a single number of 
atoms {N = 0, 1, 2), expectation values for cases where the number of fermion 
annihilation operators differs from the number of creation operators are zero. 
Using r2+(/i+) = {l + C2hj){l + cih^) and ClJ{h) = {1 + hic[){l + h2cl) we see 
that the characteristic function is of the form 

x(C,?+,ft,/i+)=Xo(?,e+)+ E X^'ii,i'^)h+h, + xT'^\i,i^)hth+h^h2 

i.j=l,2 

(46) 
and the relationship with the coefficients for ^, ^^ = 0, is 

Xo(0,0) = 1 
X2'''(0,0) = z'{dld,) 

Xf^'^(0,0) = ^'{^\^lc2^^) (47) 

Thus for one atom physical states at most six (2.2!/(2!2!) = (C^)^ + (C^)^ + (C|)2 
|40j ) c-number coefficients are required to define the characteristic function 
x(Cj C^j h, h~^) as a Grassmann function. 

The various normally ordered quantum correlation functions can be ex- 
pressed as c-number and Grassmann derivatives of the characteristic function 

G{mi,m2,n; p,l2,h) 

\d{ihty- d{ih+Y^ \d{iiY diity'^ ' 7 a(^/ll)™l9(^/^2)"^y +„„+=o 

where Tn^, 1^ = 0, 1 only, and in the case where U or mi is zero, then no differen- 
tiation takes place. 
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3.2 Distribution Function 

The characteristic function x(^, ^^, h, h~^) is related to the distribution function 
P{a, a~^ , a* , a~^* , g, g~^) via phase space integrals, which are c-number integrals 
for the field mode and Grassmann integrals for the two atomic modes. The 
formula is 



n dgtdgi I <fa+dPa 



= 1,2 

2 
X 

j=l 



expi^{gi/i+} expi{a^+} ■ P{a,a+ ,a* ,a^* ^g^g'^) •expi{^a+} expi^{/ijg+} 

i=l 

(49) 



where for the field mode we associate another pair of c- numbers a,a^, and for 
the two atomic modes i we associate another pair of g- numbers gi^gf and g = 
{(71,52}, 5^ — {.91^,52^}- As previously the characteristic function x(f, C*", /i,/i'^) 
is a function of the variables S,,^ ,hi,hf and gi,gf but not of their complex 
conjugates, whereas the distribution function P(a,a+,a*,a"'"*,5, (;"'") also de- 
pends on the complex conjugates a*,o&*. Unlike the characteristic function, 
the distribution function is non analytic. The distribution function will be a 
c- number function of a, a+, a*, a^* and a Grassmann function of gi, gt ,g2, gt ■> 
the highest order monomial being gig2g2 9i ■ The c-number integrations (Pa^ 
and (Pa are over the two complex planes a,a^ and thus Pa = daxday and 
d^a+ = daxday . The Grassmann integrals dg^ and dg^ are over the single 
variables gf and gi only, and not over gf* and 5* as well. The distribution 
function is of the positive P type for the bosonic variables a, a"*" and similar 
to the complex P type for the Grassmann variables gi,gf ■ Although pairs of 
differentials dgfdgi commute with other pairs, the convention used here is to 
write n dgfdgi = dg^dg2 dgfdgi. 

i 

The distribution function is of the form 



P{a,a'^,a*,a+*,g,g+) 
= Po(a,a^,a*,a"^*) +^P2\a,a^ ,a* ,a'^*) gig^ 

+-P2n'^^(": "'^. «*,«+*) gig2gt9t (50) 

Other possible terms do not lead to the characteristic function is given by (|46p 
when the Grassmann integration in (j49p is carried out using (I160p after expand- 
ing the Grassmann exponential factors. An ordering convention in which the 
gi are arranged in ascending order and the g^ in descending order has been 
used. The coefficients are functions of the bosonic variables a,a~^. There are 
six (2.2)!/(2!2!) = (C^)^ + (C^)^ + (C^f [40 ) c-number coefficients to define the 
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distribution function P{a,a~^,g,g~^) as a Grassmann function, the same num- 
ber of course as the characteristic. that it determines. The distribution function 
is thus an even Grassmann function of the order 2^ = 4 in the variables gi,gf, 
ff2,ff2~i ^ feature that is needed later. 

The hermiticity of the density operator leads to relationships between the 
coefficients 

Poia)* = PoC^) 

p^-\&r = P2\^) pi'\scT = pi^\^) 
p^'\&r - ^2^'(«) p^''iiir = p^''iii) 

Pr-^\ar - Pr'^\a) (51) 

where for short we write a = {a,a'^ ,a* ,a^*}. This shows that four of the 
bosonic coefficients are real and the other two P2' {a.),P2' (5) are complex 
conjugates. 

For the Jaynes-Cummings model case the quantum correlation functions are 
obtained from the characteristic function from Eq. (|48p . Applying this result by 
carrying out the differentiations on the formula (1491) relating the characteristic 
and distribution functions gives the quantum correlation functions in terms of 
phase space integrals 

G{mi,m2,n; p,h-,li) 
(fa^cfia / II dgfdgi 

x{g2t{g,Y^ {af P{a,a+ ,a\a+* ,g,g+){a+r {g+r'^{g+r- (52) 



which involves both c-number and Grassmann number phase space integrals 
with the Bose-Fermi distribution function. Note that the numbers of fermion 
annihilation and creation operators are the same. Alternative forms for the 
quantum correlation function with the distribution function as the left factor in 
the integrand and all the ai and gj placed to the right of the al and g~!f can be 
easily found. 

3.3 Existence of Distribution Function - Canonical Form 

The distribution function P{a,a^ ^a* ,a^* ^g^g'^) is not required to be unique, 
only that it generates the characteristic function x(^, ^ , /i, h'^) - which is unique 
- via the phase space integral (|1^. However, it is important to be able to show 
that a distribution function always exists. 

The existence of the positive P type distribution function can be shown by 
introducing the so-called canonical representation of the density operator, which 
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involves the fermion and boson Bargmann coherent states (see Appendix [5]) 
p = I I dg^dg / d^a+d^a Pcano7i{a,a+ ,a* ,a+* ,g,g^) A{g,g+ ,a,a^) 

dg^dg / d'^a^d^aA{g,g+,a,a+)Pcanon{a,a~^,a*,a'^*,g,g^) 

(53) 



where 



^ ^ rr(|g)5(5+*|5)Tr(|a)3(a+*|B) ^ ^ 



/« + «+* 

B 
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is a normahsed projector and 

2^ ) / / ^9^* ^9* exp( ^(g».g* + .g+*5j^ + gtgt)) exp(-- ^(aia* + a+*a+)) 

i i 

is a canonical representation function for the density operator p. The two 
forms for the density operator are the same because the normalised projec- 
tor A{g, g~^ ,a,a'^) is an even Grassmann operator and thus commutes with 
Pca7ion{ct, a'*', ct*,ce~^*,g, 9^)- By substituting for the canonical form of the den- 
sity operator ([55)) into the expression (^5)) for the characteristic function we can 
easily show that the characteristic function and distribution function are related 
as in (j49p with the distribution function given by the canonical form (j55p . For 
bosons the proof of the existence of the canonical form is given by Drummond 
and Gardiner (\15\. |18|). For fermions the proof is outlined in papers by Cahill 
and Glauber [50] and Plimak et al [21]. 

Note also that irrespective of whether or not the detailed form for P(a, a"*", a*, a"*"*, g, g^) 
is given by the canonical form (|55p . the expression ([53]) is often used to define 
the positive P type distribution function [IS], rather than the characteristic 
function expression ([49l) . If the density operator can be written as 



dg'^dg / / d^a^d^a A(g,g+,a,a^) P(a,a^, a*,a^*,g,g^) (56) 



then this defines a (non-unique) positive P distribution function. This form of 
the density operator leads to the previous relationship ([1^ between the distri- 
bution and characteristic functions. 



3.4 Probabilities and Coherences as Phase Space Integrals 

The Grassmann phase space integrals can be evaluated using the non-zero result 
/ 1 1 '^dt'^di 9^9292 9i ~ ^' giving the probabilities and coherences as bosonic 
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phase space integrals involving the six coefficients that determine the distribu- 
tion function ([50|) . 

For the two atom probability which is given by the fourth order quantum 
correlation function, we obtain the result 

Pi2^ I d^a+(faPn{a,a+,a*,a+*) (57) 

This correlation function may be zero without Po{a, a^ , a* , a^*) being zero. 
The one atom probabilities for states 1, 2 are given by 

Pi = f (fa+d^a (P2'^{a,a+,a*,a+*)- Po{a,a+,a*,a+*)) (58) 
Pa = f d^a+d^a (pl'\a,a+,a*,a+*)-Po{a,a+,a*,a+*)\ (59) 



These quantities are of course real, consistent with ([ST 
For the atomic coherences we have the results 



Pi2 = ^ d a^d aP2' (a, a^,ck*, a^*) (60) 

P21 = - f d^a+d^aP2'\a,a+,a*,a+*) (61) 



These quantities are complex conjugates, consistent with ([5T|). 
The mean number of photons is 

n= I d^a+d'a{a)Pl^'^\a,a+)[a+) (62) 

This quantity is of course real, consistent with ([?T|) . 

Note that the mean photon number result involves the fourth order expan- 
sion coefficient P4 ' (a, a"*", a*, a'*'*), and the one atom probabilities for states 
1, 2 involve the opposite second order expansion coefficients - Pj ' (a, a+, a* , a+*) 
and P2' {a,a^ ,a* ^a^*) respectively, minus the quantity Po(a,a+,a*,a"*'*) 
that determines the two atom probability. The coherences involve the second ex- 
pansion coefficients Pj' (a, a+,a*, a+*) and Pj' (a, «+, a*, a"*"*) respectively, 
multiplied by —1. Thus in their final form, only c- number integrations are 
needed to determine the quantum correlation functions. 

The normalisation integral for the combined bose-fermi distribution function 
is given by 

Trp= f f d^a+d^af[dg+dg,P{a,a+,a*,a+*,g,g+) = l (63) 



If we substitute the specific ordered form (j50p for the distribution function, we 
see that the normalisation integral gives 

d^a+d^aPl^''^\a,a+,a*,a+*) = 1 (64) 
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4 Fokker-Planck Equation 

The derivation of the Fokker-Planck equation is based on using the correspon- 
dence rules. As has been noted previously the distribution function is not unique 
and is a non-analytic function of the bosonic phase space variables. The corre- 
spondence rules are also non-unique, but the application of the standard corre- 
spondence rules (I65p - (j72p is required to lead to a distribution function which 
correctly determines the quantum correlation functions. It turns out however 
that the Fokker-Planck equation based on the standard correspondence rules 
leads to a distribution function that is not satisfactory for the Jaynes-Cummings 
model. However, the canonical form of the distribution function always exists 
and can be applied to the initial conditions. Furthermore, it turns out to lead to 
a Fokker-Planck equation that can be solved analytically so the Fokker-Planck 
equation for the canonical distribution function will now be obtained. The rele- 
vant correspondence rules will be those given in Eqs. ()81|) and the distribution 
function will be written in termsof the new variables (|78|) . 

4.1 Correspondence Rules 

We replace the Liouville-von Neumann equation for the density operator by 
a Fokker-Planck equation for the distribution function. To do this we make 
use of so-called correspondence rules, which are presented here for the general 
case of a combined system of bosons and fermions. The correspondence rules 
state what happens to the distribution function P(a, a"*", a*, a^*,^,^^) when 
the density operator p is replaced by the product of the density operator with 
an annihilation or creation operator. 
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The standard correspondence rules are: 



p 


^ 


OjP 


P{a,a+,a*,a+*,g,g+)^a,P 


(65) 


p 


=» 


POi 


P{a, a+, a\a+\g, g+) ^ L - -^"j P 


(66) 


p 


^ 


afp 


/ ^ \ 

P{a,a+,a*,a+*,g,g+)=> a+ - ^— P 

V oaij 


(67) 


p 


^ 


pal 


P(a,a+,a\a+\g,g+)^a+P 


(68) 


p 
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Cip 
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P{a,a+,a*,a+*,g,g+) ^ g,P ^ P g, 

( "d \ 
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=^ 


pn 


P{a,a+,a*,a+*,g,g+) ^ P gt ^ gt P 


(72) 


dp 

dt 


-^ 


dP{ 


a,a+,a*,a+*,g,g+) 
dt 


(73) 



The proof of these results can be obtained either starting from the expression 
(fi5|) for the characteristic function or from the expression ([55]) for the canonical 
form of the density operator. The latter approach is simpler and is based on the 
effect of the annihilation or creation operators on the Bargmann state projectors 
given in Eqs. (I188P and (|189p . Both proofs involve an integration by parts step, 
which requires the distribution function to go to zero rapidly enough on the 
bosonic phase space boundary. 

The proof of the first expressions for the fermion results do not depend 
on the distribution function P{a,a'^ ,a* ,a'^* ,g, g~^) being an even Grassmann 
function. However, the second form of the fermion results uses the feature that 
the distribution function P{a,a'^ ,a* ^a^^* , g, g^) is an even Grassmann func- 
tion. If it had been odd, then the second forms would have had all their signs 
reversed. This point becomes important when applying the correspondence rules 
to derive Fokker-Planck equation terms associated with the effects of several 
fermionic annihilation or creation operators on the density operator. The cor- 
respondence rules are carried out in succession to give the required term in the 
Fokker-Planck equation. However, multiplying or differentiating a Grassmann 
distribution function changes it between being even and being odd, and this has 
to be taken into account when applying the correspondence rules in succession 
when dealing with fermionic operators. The safest procedure is to use the cor- 
respondence rule form that does not depend on the evenness or oddness of the 
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distribution function, which is why the first forms for the fermion rules in Eqs. 
((69)) - ^ should be used. 

The non-uniqueness of the distribution function is associated with the bosonic 
phase space variables. The standard bosonic results for the correspondence rules 
in Eqs. (p5|) - (p5|) can be generalised. Since the bosonic projectors Ah{a,a~^) 
are analytic it follows that 

= I Jdg+dg J Jd^a+d^a(-^Aig,g+,a,a+)] P{a,a+,a*,a+*,g,g+) 

= I I dg+dg I I d^a+d^a (-^--^A{g,g+,a,a+)j P{a,a+ ,a* ,a+* ,g, g+) 

(74) 

and hence in the correspondence rule proof based on the canonical form of the 
density operator we see that arbitrary linear combinations 

^TTT + ^^T^ (75) 

da* daj* 

of the derivatives with respect to a* , af* may be added to each of the standard 
bosonic correspondence rules in Eqs. ([65]) - (|68)) at the applying integration by 
parts step. A similar situation occurs for the proof based on the characteristic 
function. Hence for example, the correspondence rule (|66|) can be replaced by 



p => pa^ P{a,a+,a*,a+*,g,g+)=> [a,- —— - —^:^] P = [a^- ——] P 



pdi P{a,a'^,a*,a^*,g,g~^) ^ ( a,; - -j^— p + ^ _|_^ ] P = ( "i - ^,. +, ) P 



da+* J ' V aa_, 

(76) 



(77) 



choosing A = 0, A^ = — 1 or A = 0, A^ = +1 respectively. 

There are however, further possibilities - a feature not widely commented 
upon in other work but which ultimately reflects the non-uniqueness of the 
positive P distribution for bosons. In particular, the flexibility in the corre- 
spondence rules is even greater that merely replacing a^, af, -^ or — ^ by 

these quantities plus a particular linear combination of -j^ and j_, as in (|75p 



for every term when p => di p, pa^, aj p or pdj for the creation, anihilation op- 
erators associated with a specific mode i. In fact, the linear combination used 
can be different for each a^, a,^, -^ or — ^ wherever it occurs! So for example, 

if in one term where p ^ phi we replace ( a^ — -^^ | by ( a^ — — ^ — „ i . | to 

V oaj J \ ^^i oaj J 

give I ai — _!_ J , in another term where p => pa.i we may replace I a,; — — ^ 



by [ai 
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to give ai — „, + ] ■ The reason why this is possible 



is that the only requirement is that the equation for the distribution function 
gives the correct equation for the characteristic function (or the density opera- 
tor). Additional terms of the form in (|75l) acting on the distribution function 
produce zero when in the integration by parts step they act back on either the 
exponential factor expz^jaj ^"^j and expi^{£^ja^} or the Bargmann state 

j j 

projector Ai,{a,a~^), both of which are analytic functions of the ai, af and 
hence yield zero when ^^ or _^, are applied. This flexibility is important in 

being able to convert the Fokker-Planck equation based on the standard corre- 
spondence rules into a form with a positive definite diffusion matrix (see |15] . 

For the canonical distribution function the flexibility described in the previ- 
ous paragraph is not available for the correspondence rules. They do however 
represent particular choices of the forms (j75[) . In terms of new variables 7,, 7*, 
Si, 6* which replace ai,a^ ,a*,af* via 
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(78) 



the canonical distribution function can be written in the form 



Pcanon{l-,l* ,S,S* ,g,g^) 

(4^) exp(-^(S,(S*)exp(- 



^ICiX) 



dg+*dg* exp{^{g,g* + g+*g+ + g,g+)) {g\g {j\bP\i)b \9^*) b 

(79) 



i 



The phase space integration is changed: 



(i^a+ d'^a ^ 4 



The canonical correspondence rules are 
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(80) 



(81) 
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using the results in Eqs. (|166l) . (I167p . (|169p and (|167p directly in the expression 
(ffQ]) for the distribution function. These correspondence rules can also be ob- 
tained from the standard bosonic correspondence rules by adding the additional 
terms in Eq. ([75]) . We see that 
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(82) 



(83) 

It follows that applying the modified correspondence rules in Eq. ([5^ to the 
canonical form (17^ for the distribution function gives the same result as in (I5T|) . 
since the effect of the operators involving Si, 6* is zero. Written in terms of the 
original variables ai,af ,a* ,aj* the canonical correspondence rules ()8ip were 
originally obtained by Schack and Schenzle ^S\ (see the Appendix). 
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4.2 Canonical Distribution Function - Fokker-Planck Equa- 
tion 

The Fokker-Planck equation for the canonical distribution function Pcanonij, 7*, S, S* , g, g^) 
is 

dPcanon {l , '^* , S, 6* , g, g+) 
dt 

= ^*~^ I -Q — i92)Pcano7i{l,l*,5,5*,g,g + ) + -^{gi)Pcanon{in* ^^iS* ,9,9^) 1 

+«^ ( -^™"°"('^' ^*' '^' ^*' 5' ff'^)(32')^-f + -Pcanon(7, 7*, ^, ^*, 5, 5^)(5l^)^ 1 
1 /^ , ^ * \ 

1 / ^ V \ 

+ 2*'^0 ( Pcanonil, 1* , 5 , S* , g , g^){gt)^-+ - -Pca«o«(7, 7*, f^, ^*, 5, 9^){9t)-^ 1 
-«W [-Q-^{l* Pcanon{l,l*,5,6*,g,g + )) - -^{'^ Pcanonil ,1* , S , S* , g, g+)) 

-li^(-g^^{92Kl*Pcanonil.J*,6,S*,g,g+))\ 

1 / ^ ^ \ 

+ -in {'^ Pcanon{l,J* ,S,S* , g, g+)){g} )- 



\ dgi ) 

1 fi / 1 \ 

- 2^^ ^^dijiil + ^] Pcanon{l,l*,5-,5*,g,g+)) 

+ 2*^ (5^5l)( f -Q-^ j ^canon(7, 7*, (5, (5*, ff, 3^)) 
~2*^ (((^j ^cQnon(7,7*,'5,^*,5,5^))(3^92) 



(84) 



where the bosonic correspondence rules (|8ip and the first versions of ferniionic 
correspondence rules ([55)) - ([7^ have been used. The Fokker-Planck equation 
based on the standard correspondence rules is set out in Appendix [Tni 

4.3 Coupled Distribution Function Coefficients 



Writing 7,^ for 7, 7*, ^,5 the expression ([50|) for the distribution function 
can now be substituted into the Fokker-Planck equation ()84p for the canonical 
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P+ distribution to obtain coupled equations for the six c-number coefRcients 
Po{'-f,S), P2''{'j,S) and P^ ' (7, 5) that specify the distribution function. For 
convenience we use the same terminology for these coefficients as in the general 
case and leave the canonical label understood. In the derivation the Grassmann 
differentiations of various products of Grassmann variables are first carried out 
using the results in Appendix [5] and we then equate terms involving the the ze- 
roth, second and fourth order monomials in the Grassmann variables gi, (72, gt 
and g^ to arrive at six separate coupled equations for the c-number coefficients 
Poij^S), P2"''{j,S) and P4 ' (7, ^) that specify the distribution function. 
For the zeroth order terms we have 

|Po(7, S) = -^u ( A^* _ |.^) p„(^, 5) (85) 

For the second order terms we have four equations. 



d . d \ ^i-i, 



-*^ (^9:^7* -^7 1^2^ (7,-5) 



+ \^^{l+-^)P2'\l^^)-\^^ {^* + ^)P2\l,^) (86) 



dt 



-A"ii^5) 



-*^ (a^7*-^7)^2^'(7,'5) 



-iujQP^'^ij^S) 



(87) 



Pi-\-,.s) 



-lU) 



^:l*~^l]Pr{l.5) 



+ia;oif'^(7,(5) 

+ l^nY{P2'\l,S)-P^■''i^,6)) - l^n (^|-) [p^''{^,S)^Poh,S)) 

(88) 
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+ ^zn[rP^-''{j,5)) -^^n(JP^■^\J,S)) (89) 

For the fourth order term 



9 pl2;21/ 
dt 



d ^ d X .2:21, 



1.^ / 9 A n.i/ o 1 r^ M \ nl;2, 



The corresponding coupled equations for the general distribution function are 
set out in Appendix 1101 



4.4 Initial Conditions for Uncorrelated Case 

For the case of uncorrelated initial states the density operator is a product 
'p = 'Pf'pi^ and the overall initial distribution function is just the product of the 
atomic term Pf{g,g~^) with the cavity term P{,(a, a+,a*,a+*) so that at the 
initial time the distribution function is 

Pia,a+,a*,a+*,g,g+)o = Pf{g,g+)Pbia,a+,a*,a+*) (91) 

and therefore the coefficients are given by 

Po = (d\clc2Ci\ Pb{a,a^,a*,a'^*) = 
Pl'^ = lcl^2)Pb{a,a+,a\a+*) P^'^ = - 1^%) Pb{a,a+ ,a* ,a+*) 



P 



{c\ci) Pb{a,a^ ,Q* ,a^*) Pj^'^ = {Z\cijPb{a,a^,a*,a^*) 
Pl''-^^ = Pb{a,a+,a\a+*) (92) 

where the initial distribution function for the cavity mode is 

a + a+* a* 



Pbia,a+,a*,a+*) = ( ^ j exp(-i|a-a+*n 



Pb 



2 
(93) 
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4.5 Rotating Phase Variables and Coefficients 

To proceed further it is useful to introduce rotating phase variables defined via 
the transformation 

a = /3cxp(— iwt) a* = f3* cxj){+iujt) 

a+ = ;3+ cxp{+iujt) a+* = /3+* exp(-za;t) (94) 

This then gives 

7 = —{/3 + /3^*)exp{—iujt)=jcxp{—iujt) 

7* = :^{I3* +13^) cxp{+iujt) = 7* cxp{+iujt) 

S ~ —{/3 — /3^*)exp{—iujt) = Scxp{—iujt) 

S* = ^{13* -P+)cxp{+iujt)^Tcxp{+iujt) (95) 

so that we will replace 7, 7*, (5, 6* by 7, 7*, (5, (5 , or 7, (5 for short. The canonical 
distribution function coefficients will be of the form P(j,S), where we now 
change the labelling of the coefficients to P so as to reflect this variable change. 
This transformation enables the elimination of the cavity field term from the 
equations for the coefficients. Explicit time dependences will be left understood. 
For any P(7, 6) 

|^(7,'^)=»^(^7m^)-^rm^)) +1^7,^) (96) 

where since each canonical distribution coefficient is of the form exp{^6 S ) F{'j,'-f ) 
it follows that the derivative terms involving i5, S can be ignored because 

Note also that the original 7, 7*, 6, 6* were time independent. On the right side 
of each coefficient equation we have a term 

This cancels out the equivalent term on the left side so we now have the following 
results. 

For the zeroth order terms we have 

^^0(7,^) - (99) 
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For the second order terms we have four equations. 



-iVt ( 7* + ^ j -^2 '^(7 J) exp(+ia;t) 



-2WoP2'^(7,(5) 



,^7 



2 

izf} exp(-zc.t) ('^') (P2'''(7 J) - ^0(7,^) 



+2WoP2^'^(7,(5) 



~ izr! cxp(+zc.t) (^^^ (P2^'(7,^) - ^0(7,^) 



|i?^^(7J) 



1 / ^ ~i 2 ~ ~ 
+ 2*^(7 ■P2' (7, <5) cxp(+ia;t) 

- 2*^ (7-^2^'^ (7 J) cxp(-ia;t)j 



For the fourth order term 

d J512;21/~ Tn 



dt 



1 .„ / 9 \ ~2:1,~ 7, , . ,x ^ cf ^ \ nl:2 



(100) 



(101) 



(102) 



(103) 



-i^ [^)P2 (7,'^)exp(-icjt)- -ir^ ( — j P2' (7, f^) cxp(+iwt) 

(104) 
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Also, equations with nq^expUcit time dependence can be obtained via tjie 
change of coefficients from P to 5 by incorporating exp(±iwi) in the Pj ' (Tj '^)) 
P2 ' (7, S) coefficients, and we can also factor out the explicit dependence on 6 
via the overall factor exp{—6 5 ) = exp(— 5^*). We also note that the canonical 
form (|75|) of the distribution function (and hence all the coefficients) contains 
the factor exp(— 77*) = exp(— 77*), so it is useful to eliminate this factor via 
the change from P(7, S) to S{j). The new coefficients are: 

P2'^{l,'S) = Sl'^il) exp(-77*) exp(-'^^ ) 

P2''^{l,'S) = ^2'^(7) exp(-77*) exp(-'^^ ) 

P2' (7,(5) — S'2' (7) exp(— 77*) exp(— (5(5 )exp(— iwi) 

P2'^ilJ) = 52'^(7) exp(-77*) exp(-(5(5 )exp(+ia;i) 

Pq{i,S) = 5o(7) exp(-77*) exp(-(5(5 ) 

Pf^2i(7,^) = ^f ^'^(7) exp(-77*) exp(-^?) (105) 

The separate cavity and atomic transition frequencies are then incorporated 
into the detuning 

A^uo-Lo (106) 

Also, S now plays no further role in the dynamics and the S only depend on 
7,7* (or 7 for short). 

With this substitution the equations for the coefficients are as follows: 

For the zeroth order term 

^^0(7) = (107) 

The four second second order equations are. 
d 



^^{sl'\^)-So{i)) 



dt 



l^n (^) sl'\^) - Un f A ) s},'\^) (108) 



^ cl:2/~N 



Ua^ [S'^\^) - 5o(7)) + ]^^ (^) (5^(7) - ^0(7)) (109) 






\^nr{s',^\l)~S,{^)) - \^^ (^) (^2^'(7)-5o(7)) (HO) 



2 
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1(5^(7) -5o(7)) 
= +l^n[rsl■''{l)) -l^n{JSl■'\J)) (111) 

where we have substracted the zero quantity ■4fSo{'y) from each side of the first 
and fourth equation. From subsection 13.41 we see that the quantities S'2*(7) — 
50(7) determine the one atom probabilities. However, from ()107p and the initial 
conditions (l92|) we see that <S'o(7) will be zero for the one atom Jaynes-Cummings 
model, so we can ignore 5*0(7) henceforth. 
The fourth order equation is 



9 pl2;21 

dt 



S4 (7) 



1 .^ /^ ~ 9 \ ~2:l,~s 1 -r^ /^ ~* , ^ \ ol:2 



Analogous equations for the standard distribution function are set out in Ap- 
pendix [TUl 



5 Solution to Fokker-Planck Equation 

The approach used by Stenholm [13] can be adapted to provide an analytical so- 
lution for the canonical distribution function for the one atom Jaynes-Cummings 
model for any initial conditions. The equations (|108p - (lllip for the 82'' (7) can 
be solved via the substitution 

^2^''(7) = **(r)^,(7) (113) 

where the ^^ (7*) are functions of the 7* and the ^4(7) are functions of the 7, 
and where ^4(7) satisfy the coupled equations 

-^1(7) = -iA*l(7)--^^! (^— j*2(7) 

-^2(7) = -iiA^2(7)-^*^7*i(7) (114) 

This ansatz is consistent with the original equations (|108p - (|llip for the S^''' (7) . 
As indicated previously we have set S'o(7) = for the one atom case. 

Differentiating the second equation and substituting from the first gives 

^*^(7) + Ja^*2(7) = -i^'7(^)*.(7) 
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where the substitution 

s = Ig 7 7 = exp s (116) 

has been made. 

A solution of the equation (IllSp can be obtained using separation of the 
variables 

^2{l) = T{t)Kis) (117) 

whence we find that 

where since the left side is a function of t and the right side is a function of s 
the quantity A must be a constant. 

The solution of these two equations is straightforward. We have 



K{.s) = Cexp(-js) 



4A 

(4A/ n'^) 



= C(7)^""'"^ (119) 

where C is a constant, and 

T{t)= Acqs (^/xTAylt) +B sin (^/xTA^t) (120) 



with A and B also constant. 

Since we require the overall distribution function to be a non-singular single- 
valued function of the phase space variables we see from Eq. (I119|) that there is 
a restriction on A such that 

4A 

— =n (n = 0,l,2,..) (121) 

where n is an integer. 

Combining the variables to eliminate A and absorbing C into the other con- 
stants we see that a solution for ^2(7) for a particular integer n is 

^2(7) = 7" f^n cos -Lu„t + S„ sin -LoM (122) 

where 



uJn = \Jn n^ + A2 (123) 

is the frequency associated with population and coherence oscillations in the one 
atom Jaynes-Cummings model. The corresponding solution for ^1(7) is then 
obtained from Ea. (jll4l) and thus 

.(„_i) f__^l_ , { ujnBn + iAAn] , ^ 1 ^ f -a;„A„ + iA_B„ 



$1(7) = iY [ cos -uj„t< ^+sin-a;„i 



2 " 1 n I 2 " 1 n 

(124) 
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However, a solution with n = leads to a singular 7" behaviour, so it follows 
that n is restricted to the positive integers. Also, as the ansatz equations are 
linear the general solution is a sum of terms with differing n so that we finally 
have the solution in the form 

*i(7) = *1,J* '[[ ^ |cos-..„t+| |sm-c.„t 



*2 (7) = X! '^'" [ ^" "^^^ -UJnt + Bn sin -UJnt j 

The constants An, Bn are chosen to fit the initial conditions. 

We can now express the original distribution function coefficients in terms 
of these function using Eqs. (|113p and (|f05p . Expressions from which S^ ' (7) 
could be obtained, but these are of little interest. The results can be written in 
terms of the original phase variables by substituting 

J ^ -{a + a'^*)exp{+iujt) 2 = -{a - a^*) exp{+iujt) (126) 



from (|94l) and (|95|) into the above results. 

The analogous solution for the general distribution function is set out in 
Appendix [101 It turns out that the Fokker-Planck equation based on the stan- 
dard correspondence rules leads to equations for the coefficients that can also be 
solved by a similar ansatz. However, the solutions lead to a distribution func- 
tion that diverges on the phase space boundary and in general diverge at large 
i, with dependences on hyperbolic functions of ^V^^t in the case of zero de- 
tuning. This then throws the original derivation of the standard Fokker-Planck 
equation into doubt because the integration by parts step fails. Other cases 
where this occurs have been studied by Gilchrist et al [41]. 

Finally, as we will see in subsection 15. 1[ the solutions based on the canonical 
distribution function (|105p agree with the standard quantum optics result, and 
in particular the quantities A^, B^ can be chosen so that the initial conditions 
are the same as for the canonical distribution function determined from the 
standard quantum optics solution. 



5.1 Comparison with Standard Quantum Optics Result 



As a comparison, we now calculate the canonical distribution function (|55p 
as determined from the state vector given in Eq. ^ and ([5]) obtained from 
standard quantum optics methods. 



(125) 
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The density operator is 

n,rn 

+ E B^^2it)B^i2it) exp(-z^^^T^c.i) 4 |0) (0| ?2 j^T^T) (;^r^T| 

n,Tn 

+ E 5;i3T2(i)5™i(i) exp(-z7r^^u;t) ?t |o) (0| ?i [^T^T) (m| (127) 
where the new amphtudes are given by 

S„i(i) = [cos(-ujrit){A„i[0)} + ism(-uj„t) ' 



5;r:T2W == ( cos(-a;„t)|A;j^2(0)| -«sin(-w„t)| 

(128) 
noting that 

exp(— -iAf) exp(— ?(naj — -wo)^) = exp(— i(n — -)a;i) 
exp(+-iAi)exp(— z(n — Icj + -a;o)t) = cxp(— i(n — -)wi) (129) 

Then with A = ^{(^ + Q^^*); ^ — ^{o/- ~ ct~^*) and 7 — 'j exp{+iujt) , S = 
6 exp{+iLdt) , as in Eqs. (l78l) and (|95|) we have on substituting into Eq. (|55l) and 
using expressions for the bosonic Bargmann states from Appendix [9l 



B\ -{a + a+*) 






j2B„,it)B*ut)c\\o)mi (^*)"(^)''' 



' \/m} 



■Y.B—At)B^,,{t)c\ 



n\ ymi 

*\n — 1 /~\m— 1 



C2 



(7 )""^ (7)^ 



VF^WWirr^-^V- 



and 



+ ^i?„i(t)i31^2W cl |0) (0|C2 (7T (TJ:i^ exp(+»^0 

n,m ^ V V / 

+ T.B—i2(t)B*Ut) c\ |0) (0|ci ^lZl_ ^exp(-zcut) (130) 

^ V(n-l)!Vm! 



j^ * 

exp(— -(aa* + a+*a^) = exp(— 77*) exp(— (5(5 ) (131) 
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Also using expressions for the fermionic Bargmann states from Appendix [3] 

B{g\ (ct|0)(0|cfe) \g)s 



^9t9p StiSjk 



(132) 



so that 



ij 



B{g\ B (-(« + «+*) 



-(a + «+*)) \g} 



j2Bni{t)B:^,it) 



*\n (~\m 



9i9i 



irr (7) 



' A/rr?' 



\("-i) r^^^C"'-!) 






v/(n-l)!v/(™-l)! 
+ ^i?„i(i)i?l^,(i) glgt* ^Vi /T ,„exp(+»c.O 

+ E^^2W5:a(i) ff25r ^V ... T . expHc^t) (133) 



\/{n — 1)! vto! 

The Grassmann phase space integrations can be carried out using results 
from Appendix [8] and give 

/ / dg+*dg* exp( ^(3^5* + gt*gt + g^gt)) algt* 

i 

= gwt {5k25n) + gigt {-5k25ii) + g2gt {-5ki5i2) + 525^ {5ki5ii) 
+gig2gtgt i^k25i2 + SkiSii) (134) 

so combining the results we find that the canonical distribution function is 
Pcanon{a, a'^ , a* , a'^* , g , g'^) 



V47rV 



exp(— 77*) exp(— (5(5 ) 



[J2 Bni {t)B*^, it) {g2gt + gig2gt gt } 



*\n (~\m 



irr (7) 



-J2B^^,^it)B^^{t) {gigt+gig2gtgt} 



71! y/ml 



V(n-l)! V(m-l)! 



^/n\ \/(m — 1 
.J2B—,{t)B*^,{t) {-g,g+} ^J/ 7 exp(-»c^t)] (135) 

n.TTi V V y V 
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From this result we can identify the coefficients 



(T) 



(n-l) 



n,m 

, ^^2 , exp(-7r) exp(-^r)^i?— 2(t)B,;;,i(i) 



(7) 



(m-l) 



-P2'^(7,'5) = -exp{-iujt)( 



V(n-l)! V(m-f)! 



i'2^'(7,^) 



exp(+iwt) 



Utt^ 



exp(-77*) exp(-<55*)^i?„i(i)i?^2W 

(7*)" (7)" 



■\/(n — f )! Vm! 



^'^'(X-^) = f ^ j exp(-7r) exp(-<5<5 )5] i3,aWi?™i W ^= 



n\ y/{m — 1)! 
(136) 



and 



Po(7,<J) = 



Pf^2i(7,<5) 



1 

47r' 



2 , exp(— 77*) exp(— (5 5 ) 



x[^i3„i(i)i?;„i(i) 



(7*)" (7)'^ 



J2B—uit)B^uit) 



(137) 

(;y*)("-i) (;y)(™-i) 



If we write 



$1(7) 
^2(7) 



I 



-]yB. 



it) 



(7) 



(m-l) 



2nJ ^ m-12V ^ ^(^_1)! 

(7)" 



i E^^iW 



(138) 



then 



P2'^i(7,^) = exp(-77*)exp(-JO'i'i(7*)$i(7) 
P2^'^(7,I) = exp(-77*) exp(-jr)$J(7*)$2(7) exp(-iwi) 

^2 (7i'^) = exp(— 77*) exp(— (5(5 ) $2(7*) *^i(7) exp{+iujt) 
P^^ilTS) = exp(-77*)expM?)$*(7*)$2(7) 

To have agreement with the previous results in Eq. (IIOSP so that 

$1(7) = ^1(7) 
^2(7) = ^2(7) 



(139) 



(140) 
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where ^1(7) and ^2(7) are as in Eqs. (|125p . we require 



n 



Bn — 





4^1^2(0)}*^ 


1 


rT2(0)r 






n-l)\ 

m+AA- 


1 


^ 

~){A„i(0)}*- 

V 


1 


J 

^—2(0; 




i-iy. 

1 


)\'^) \ 


OJn 




(141) 



The last two equations give explicit expressions for An and B„. Substituting 
these expressions into the left side of the first two equations gives the right hand 
sides, showing that the four equations are consistent. Hence we see that for any 
initial conditions for the one atom Jaynes-Cummings model, the solution given 
by the Grassmann phase space approach is the same as that from the standard 
quantum optics treatment. 



5.2 Application of Results 

As an illustration of how to apply the above results for the canonical distribution 
function we consider the case where the atom is initially in the lower state and 
the field is in a coherent state of amplitude 77. In ths case we have from subsection 
7? 



p2;i _ n 



3i;2 
2 

32;2 



pr = 



P^'^ =Pbia,a+) 



(142) 



with 



n(a,a^ 



47r2 



exp(- 



\a — a 



+ *|2 



-)exp 



-{a + a+*) — T] 



Hence 

P2'^{l,'S) = ^ exp(-'^^ ) exp (- I7 - r?|^) = cxp(-'^^ ) exp(-77*) S^''^{j) 



i-f 



1 



'S'2' (7) = TT°^P(^(7-'7)(7*-?7*))exp(+77*) 



47r2 
1 



2tj. cxp(- 2 W*) exp(-7*77) 



— exp(--?7?7*) exp(-777*) 

Ztt Z 



(143) 
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At i = we have 



*i(7) = ,^^(n^i)(J -n^^ + ^AA. ,, ^ 



n=l 



_ °° _ 11 

^2(7) - ^7"(^n) = ^exp(--r77?*)exp(-7r/*) (144) 

n— 1 

SO that if we choose 



1 / 1 *A-v*y 



27r "^ 2" ^ n! 
Bn = ^A„ (145) 

the solutions for ^1(7) and ^2(7) determines the time dependent distribution 
function. 



6 Conclusion 

We have shown that a phase space approach using Grassmann variables to de- 
scribe the atomic system and c-number variables to describe the cavity mode can 
be used to treat the Jaynes-Cummings model and to obtain the same results for 
treating phenomena such as pure Rabi oscillations and collapse, revival effects 
as those from standard quantum optics methods. The Liouville-von Neumann 
equation for the density operator was converted into a Fokker-Planck equation 
for the canonical positive P distribution function using the correspondence rules 
associated with this choice of distribution function. The distribution function 
is a Grassmann function involving Grassmann phase space variables 51 , g^ and 
92i 92 fo'^ the two fermionic modes associated with the two atomic states, with 
six c-number functions of the bosonic phase space variables a, a+ associated 
with the cavity mode being involved as coefhcients in specifying the distribution 
function. In the context of a general mixed state where there may be zero, one 
or two atoms present, expressions for the probabilities of finding one atom in 
one of the two atomic states, one atom in both atomic states and no atom in 
either atomic state were obtained as bosonic phase space integrals involving the 
six bosonic coefficients. Coupled equations for the six bosonic coefficients for 
the canonical distribution function were obtained from the Fokker-Planck equa- 
tion These equations were solved for the one atom Jaynes-Cummings model 
using an ansatz similar to that applied by Stenholm |14| in an earlier Bargmann 
state treatment of the Jaynes-Cummings model, and the results shown to be 
equivalent to the standard quantum optics treatment based on state vectors and 
coupled amplitude equations. 

Positive P distribution functions have the feature of being non-unique, with 
different correspondence rules applying in the derivation of the specific Fokker- 
Planck equation. In this application we also found that applying the standard 
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correspondence rules (rather than those for the canonical positive P case) leads 
to a Fokker-Planck equation where the solution for the coupled equations for the 
bosonic coefficients via a similar ansatz was quite unsuitable. Not only did the 
solutions diverge for large time t, but the distribution function diverged for large 
phase space variables a, a^, thereby throwing into question the derivation of 
the Fokker-Planck equation. The standard treatment requires the distribution 
function to vanish on the phase space boundary. As the correspondence rules for 
the canonical positive P distribution do not require this feature, it is suggested 
that Fokker-Planck equations based on the canonical positive P distribution may 
be more reliable. Furthermore, in terms of matching a general solution of the 
Fokker-Planck equation to the initial conditions, the use of the canonical form 
of the distribution function is easiest since initial conditions are usually specified 
via the initial density operator from which the canonical distribution function is 
directly determined. However, more general Fokker-Planck equations involving 
derivatives higher than second order may occur using the canonical distribution 
function |28| . so that no replacement by Langevin stochastic equations is then 
possible. In fact even if the Fokker-Planck equation is only second order, the 
diffusion matrix may not be positive definite. 

The successful treatment of this classic quantum optics system based on 
phase space methods using Grassmann variables represents an important step 
in applying such methods to treat more complex problems involving fermionic 
systems. 
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8 Appendix 1 - Grassmann Numbers and Cal- 
culus 

Since Grassmann variables and their calculus may be unfamiliar to many physi- 
cists a short summary may be desirable. More extensive accounts of the prop- 
erties of Grassmann variables and their calculus are given in Refs. ^01 [TO] . 

8.1 1.1 Grassmann Algebra 

Grassmann variables satisfy the following anti- commutation rules with each 
other and with fermion annihilation, creation operators 

9t9j = -gjm 9iCj = -Cjgi gic] = -c^g^ 
{91:9]} = {9i,cj} = {gi,cl} = (146) 

and Grassmann variables and fermion operators commute with c-numbers and 
boson operators. A key feature of Grassmann variables that immediately follows 
is that their square and hence all higher powers are zero. 

fff = <?f = ... = (147) 

Grassmann variables also have no inverses, and hence division is undefined. 

Grassm,ann functions involving linear combinations of products of several 
Grassmann variables with c-numbers coefficients can be defined in an obvious 
way and are of the form 

/(/ii,/i2,..,/i„) = /o + S /t ^j + E /y ^i'ij + E fijkhihjhk + .. 

i i<j i<j<k 

+ fl23..nhih2..hn (148) 

where /o, fi, fij, .., /i2..n are c-numbers. Functions may be even or odd depend- 
ing on whether all their terms contain an even or odd number of Grassmann 
variables. A general function is the sum of an even and an odd function. Two 
even or two odd Grassmann functions commute with each other, whilst an even 
and an odd function anti-commute. 

Grassmann functions have a linearity feature in that each Grassmann vari- 
able can at most appear in a linear form due to (|147|) . Thus 

/(/ll,/l2,..,/l„) 

= at{hi,..,hi-i,hi+i,..hn) + biihi,..,h^^i,hi+i,..hn)hi (149) 

= ai{hi, .., hi-i, /li+i, --hn) + hi Ci{hi, .., /ii_i, ft-i+i, ../in) 

where at , bi and Ci are Grassmann functions that are independent of hi . In gen- 
eral hi and Ci are not the same. Linearity has the effect of truncating functions, 
thus the Grassmann function exp g is just equal to I + g. 

Complex conjugation can also be defined with the property 

{g.gjr = g^g: (150) 
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Apart from the conjugation, the anti-commuting multiphcation rules and the 
lack of an inverse, Grassmann variables satisfy the normal algebraic rules (as- 
sociative laws of addition, multiplication etc) . 



8.2 1.2 Grassmann Differentiation 

The basic rules for Grassmann differentiation are 

|-5. = 5.|-=^.. (151) 

where both left and right differentiation are defined. 

Differentiation proceeds via moving the variable to be differentiated to the 
left or right of all the other Grassmann variables and then applying the above 
rules. Thus for the Grassmann function in (J149I) we have 

^Ti-f{hi,h2,--,hn) = Ci{hi,..,hi_i,hi+i,..hn) 
ohi 

f{hi,h2,--,hn)j—- = bi{hi,..,hi-i,hi+i,..hn) (152) 

Note that in general, left and right differentiation do not give the same result. 
Multiple differentiation is carried out in the order of the derivatives. Thus 
for left differentiation 

t ^ -^ (^ \ 

with equivalent results for right differentiation. Mixed left and right differenti- 
ation also occurs and it turns out the different orders for carrying this out give 
the same result. Also left and right differentiation are related for even and odd 
functions. 



jrMg) = {-l)fE{g)j^ (154) 

l^Joig) = (+l)/o(.)| (155) 

Product rules for differentiation can be derived. These depend on whether 
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the factors are even or odd Grassmann functions 



ly^'-'^ -' 


- ^l'^>f'^f'^al/'> 


%''?« = 


- ^lj?^''?i,'^> 


'"^'l - 


- «^.%;) + <^^al)^f 


<"'>l - 


- M''a,?-'J'ay? 



(156) 

Thus the product rule is different in general from that in ordinary calculus. For 
Grassmann functions that are neither even nor odd the derivative of a product 
can be obtained from (|156p after writing the function as the sum of its even and 
odd components. 

8.3 1.3 Grassmann Integration 

The basic rules for Grassmann integration are 



dgil = Idgi = 

dgigj = Sij. gjdgi=~Sij (157) 



where both left and right integration are defined. The different results for left 
and right integration are due to the differentials thems being anti-commuting 
Grassmann variables. In the present paper only left integration will be used. 

The same outcome from -^gj = Si,j, -^1 — and / dgt gj — 5ij, / dgil = 

provokes the comment that differentiation and integration are the same, but 
this is not really the case in view of the different results for right integration 
and differentiation. 

Integration proceeds via moving the variable to be integrated to the left or 
right of all the other Grassmann variables and then applying the above rules. 
Thus for the Grassmann function in (J149I) we have 



dhif{hi,h2,..,hn) = Ci{hi,..,h,_i,hi+i,..hn) 
f{hi,h2, .., hn) dhi /(/ii, /i2, .., K) = -h{hi, .., /ii_i, /i^+i, ../i„) (158) 

Multiple integration is carried out in the order of the differentials, for exam- 
ple in left integration 



dhidhj f{hi,h2,..,hn) ^ dhi dhj f{hi,h2, ■.,hn) ] (159) 
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with equivalent results for right integration. 

An important results is for the complete Grassmann integral 



.Jdh,dh,..dh,.h,....h,h,^l (160) 



8.4 1.4 Grassmann States and Grassmann Operators 

For fermion systems, vectors can be defined in a generalised form of Hilbert 
space, which involve linear combinations of basis vectors such as Eq.(??) but 
now with Grassmann numbers as the coefficients. Such Grassmann vectors 
are not taken to represent physical states (where the coefficients must be c- 
numbers) even if a fixed number of fermions are involved, but they have uses in 
the mathematical manipulations. The fermion coherent states are Grassmann 
vectors. 

Similarly, we may introduce generalised operators in this new Hilbert space 
by taking linear combinations of the products of the fermion operators with 
Grassmann numbers as the coefficients, such as in Eq. pS)) . Such Grassmann 
operators do not represent physical quantities (which must involve c-numbers 
as coefficients) or symmetry operations, but again are useful mathematically. 

Many of the results for Grassmann functions also apply for Grassmann op- 
erators and states, in which some of the Grassmann variables are replaced by 
fermion annihilation or creation operators. The anti-commuting feature of the 
fermion operators with Grassmann variables enables the same proofs to be made. 

In addition to the rules (jl46p for Grassmann numbers, it is necessary to state 
the basic rules for multiplying the vacuum state with Grassmann number, since 
the Fock states involve products of creation operators acting on the vacuum 
state |0), . The rule is that the Grassmann numbers commute with |0) or (0|. 

g\0) = \0)g .g(0|-(0|g (161) 

Note that as a consequence 

g\mi;m2]n) = (-l)™i+'"" |toi; ma; n) 5 

(mi;m2;n|5 - (-l)™i+"'^g (rm; ma; n| (162) 

so clearly for N fermion states, a Grassmann number anti-commutes with the 
Fock states for fermion systems with an odd number of fermions, and commutes 
if the number of fermions is even. The bra and ket vectors for fermion states 
may be classified as even or odd vectors depending on whether they only contain 
terms with even or odd numbers of fermions. Thus (|0; 0; ri) + |1; 1; n)) would be 
an even vector, whilst (|0; l;n) + |l;0;n)) would be odd. These concepts may 
be extended to include Grassmann vectors as well as ordinary state vectors. 
In this case as well as the vacuum state |0) or (0|, all terms in even (odd) 
vectors contain an even (odd) number of fermion creation operators, annihilation 
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operators and Grassmann numbers, and consequently commute (anti-commute) 
with any Grassmann number. Thus the Grassmann vectors (1 + cj/ii) |0;0;n) = 
\0;0;n) — hi |l;0;n) or (0;0;n| {1 + hfci) = (0; 0; n|~(l; 0; n| /i^ are even vectors. 
Some operators such as those that represent physical quantities contain an 
even number of ferniion creation and annihilation operators, usually the same 
number of each (sec the Hamiltonian H in (j26p for example, which is the sum 
of terms each containing the same numbers of annihilation and creation opera- 
tors). These operators are called even operators. Consequently even operators 
will commute with a Grassmann number. Other operators of interest such as 
the ferniion cretaion and annihilation operators themselves involve odd numbers 
of creation and annihilation operators. These operators are called odd opera- 
tors. Operators either commute or anti-commute with a Grassmann number 
depending on whether they are even or odd. Arbitary operators can always be 
expressed as the sum of an even operator and an odd operator and a Grassmann 
number will commute with the even component and anti-commute with the odd 
component. These concepts may be extended to include Grassmann operators 
as well as ordinary quantum operators. In this case all terms in even (odd) oper- 
ators contain an even (odd) number of ferniion creation operators, annihilation 
operators and Grassmann numbers, and consequently commute (anti-commute) 
with a Grassmann number. Thus the operators tlt{h'^), QJ (h) in (|45|) are even 
operators. Operators for which there are no Grassmann variables involved are 
just a special case where evenness or oddness only depends on the number of 
fermion creation, annihilation operators. 
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9 Appendix 2 - Bargmann Coherent States 

For both bosonic and fermionic systems we can define the coherent states, which 
are defined by the effect of unitary displacement operators on the vacuum state. 
These states are parameterised via c-number variables ai,a* for the bosonic 
modes and via Grassmann variables gi,g* for the fermionic modes. The coher- 
ent states are normalised to unity and are eigenstates of the boson or fermion 
annihilation operators with eigenvalues a^ or gi respectively. A full description 
of the coherent states is given in [42], [20] for the bosonic and fermionic cases. 
As in these papers, we will treat the general case of multi-mode systems. For 
our purposes it will be convenient to use a related set of un-normalised states 
called the Bargmann states, which have the property of only depending on a^ 
or gi, and not on the complex conjugates. 



9.1 2.1 Bosons 

For bosons the Bargmann states are defined as 

\a)g = exp r£alaA |0) =l[exp (ala,) |0) (163) 

which only depend on the c-number variables a = {ai, a2, .., a^, .., ««}, and not 
on the complex conjugates a* = {a|, a2, ■■, a*, .., a* }. These are related to the 
coherent states \a,a*) via 



\a,a 



exp(--a* ■a)\a)g 


(164) 


exp(-la*.«)n>'^h), 


(165) 



where a* ■ (3 — Y^ cutPi- The second expression is the well-known expansion of 

the coherent state in terms of Fock states. This c-number expansion represents 
a Poisson distribution of number states with a mean boson number (n^) = |aip 
and a variance (Anf) = Ijii). The coherent state is not a physical state except 
in the case of photons. 

As in the case of coherent states, the Bargmann states are eigenstates of the 
annihilation operator 

ai\a)g = ai\a)g (166) 

B{a\a\ = B{a\a*- (167) 

They satisfy normalisation and orthogonality conditions 

B(a|/3)B=exp{a*-/3} (168) 



Hence Bargmann states are unnormalised versions of the coherent states. 

The operation of a creation operator on a Bargmann ket vector or an annihi- 
lation operator on a Bargmann bra vector can be written in terms of derivatives 
of these vectors. 



«1 \^)b = [^] W)b (169) 



d_ 
da 

^) ("Ib (170) 

there being no distinction between left and right differentiation. 

9.2 2.2 Fermions 

For fermions we define the Bargmann states via 

\g)^ = exp ij^ c\gA |0) = X{{1 + ig^ |0) = X{{\Q,) - g, |1.)) (171) 

\i— 1 / i i 

which only depend on n Grassmann numbers g = {gi,g2, --i diT-j Qn}, and not 
on the complex conjugates g* = {gl,g2, --,9*, ■■^dn}- These are related to the 
fermion coherent states \g,g*) in a similar way as in the boson case. 



\g,g*) = exp(--g*.g)|g)5 (172) 



exp(-ig*-g)n(|0':)-5dl.))- (173) 



where g* • h = ^ 9* hi. This differs from the corresponding bosonic expansion 

because it only involves a superposition of a zero fermion state with a one 
fermion state. This is to be expected from the Pauli exclusion principle since 
any mode can only be occupied by at most one fermion. The fermion coherent 
state is of course unphysical as it involves Grassmann numbers as expansion 
coefficients. The Bargmann states are employed by Plimak et al. [2T] rather 
than the fermion coherent states as in Cahill and Glauber PD|. Bargmann states 
are even Grassmann vectors. 

As for the coherent states, the fermion Bargmann states are eigenstates of 
the annihilation operator 

cdg)s == .9*1.9)5 (174) 

B{g\cl = B{g\9*. (175) 

In the fermion case we can also find eigenstates of the creation operator, see 
[20] . The Bargmann states satisfy normalisation and orthogonality conditions 

B {g\h)g = exp{g* • h}. (176) 
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Note the similarity of these results to those for the boson states. 

The operation of a creation operator on a Bargmann ket vector or an annihi- 
lation operator on a Bargmann bra vector can be written in terms of derivatives 
of these vectors. In the fermion case 



d_ 



^1\9)b = {-7r-]\9)B = \9)B{+—] (177) 



!-] = ( 1. 

99* Vdg* 



{9\b c^ = {9\b I -— I = I +— I {9\b (178) 



where we note that both left and right derivatives forms apply for each of cj \g)g 
and {g\g Ci. This is analogous to the feature that in (I174p and (|175p the eigen- 
values can be placed on either side of the bra or ket Bargmann vector. 

The Bargmann states can be used for a representation of quantum operators. 
For the fermion operator Qf we introduce completeness relationships for two 
sets of fermion Fock states \vi^V2, ■■■,Vn)-, \^iT^2T-iCn) j where the occupation 
numbers for the various single particle states Vi, ^j are 0, 1 only. We have 

{9\B^f\h)B^ Y^ Y^ (5|l^l,1^2,-,I^n)B(^l:i'2,-,l^«|%|Cl,?2.->C„) (Ci,C2>-'C„|/i)b 

Then from the expression for the Bargmann states 

(Ci, ^2, -, uh)B = n (^j 1 11(10^) + ii*) ^^) ^ n(^«.o + -^e^i^*) = ne^^)*' 

j i i i 

giving the result 

{g\gClf\h)g= J2 E ^fiM■A^}){9LJ'i9:nJ'■■{9*raJHhn^J".■{h^2)^Hh^,)^^ 

(179) 
where nf{{iy};{^}) = {vi,V2, ■.,Vn\^f\^i,C2, -^D- 

9.3 2.3 Projectors 

The Bargmann states can be used to define normalised projectors as 

where the trace of the projectors \a)g {f3\g and \g) ^ {h\g have a simple form 

Trf{\g)B{h\B) = cxp{g ■ h*} =b {-h\g)s =b {h\ - g)^ (182) 

Tn{\a)s{l3\B) = cxp{a • /3*} =s (/3|a)B (183) 
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The normalised projeetors have the property that their trace is unity. 

rrfcAb(a,/3*) = l Trilf{g,h*) = l. (184) 

9.4 2.4 Completeness 

One of the important features of both the boson and fermion coherent states is 
that they satisfy completeness relationships. In terms of Bargmann states the 
completeness relationships are 

J (fgexp{-g* -g) \g)s{g\B = 1 

/ (Pa.cxp{—cy.* ■ a) Iq;)^ (aj^ = ttI (185) 

for fermions and bosons respectively. Here d^g = J^ dg*dgi and d^ct = J^ dui^daiy 

i i 

9.5 2.5 Trace Properties 

Cyclic properties of Grassmann operators can be established. For a pair of even 
Grassmann operators 

Tr{Clf{h)Af{k)) = Tr{Af {k)nf i-h)) = Tr{Af {-k)Clf (h)) (186) 

For the product of two odd Grassmann operators. 

Tr{h^{h)K^{k)) = Tr{k']{k)h']{-h)) = Tr{k^ {-k)h^ {h)) (187) 

9.6 2.6 Operator Identities 

The effect of bosonic annihilation and creation operators on the Bargmann state 
projectors involved in the canonical form for the density operator are 

a| Ab(Q:,Q:+) = (- ^ al j Ab{a,a'^) 



( — 
\da1 



Ab(a, a+)ai = [ ^^-^ + ai] Ab{a,a^) (188) 
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The effect of fermionic annihilation and creation operators on the Bargmann 
state projectors involved in the canonical form for the density operator are 

CiA/(g,g+) = giAf{g,g+)=Af{g,g+)gi 
A/(5,.9+)cl = Af{g,g+)g+ ^g+Af{g,g+) 

^f{9,9^)c^ = A/(.9,.9+) I -^-.g,J = I +^-.g,J A/(.g,.g+) 

(189) 

the second form of the result following from Af{g, g^) being an even Grassmann 
operator. 
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10 Appendix 5 - Results for Standard Distribu- 
tion Function 

10.1 5.1 Fokker-Planck Equation 

The Fokker-Planck equation for the general distribution function P{a, a^ , a* , a^* , g, g~^) 
is 

dP{a,a+,g,g+) 
dt 

Ea ( ~^ ^ \ 

= -«-^ ( — {32^(a,a+,5,5+)} + — {5iP(a,a+,g,g+)} j 

+,^ l{P(^a^^+ ^g^g+)g+}A^ ^ {P{^^^a+ ,g,g+)gt}^\ 

"2*^0 ( Q—{92Pia,a+,g,g+)} - —{giP{a,a+,g,g+)} j 

1 / V V \ 

+ -iujo ( {P{a,a+,g,g+)gt}-—^ - {P{a,a+ ,g,g+) g^}—^: 1 



g^{aP{a,a ,g,g )} + — 



-luj I -^{aP(a,a+,5,.g+)} + ——{a+P{a,a+,g,g+)} 



1 (^ ^ \ 

~ 2* " I a~^^2a^ ^(a, a+, .9, .g+)} + -^{gia P{a, a+ , g, g+)} 1 

if V V \ 

+ -in {Pia,a+,g,g+)g+ a}—— + {P{a,a+,g,g+)g^a+}--r 

-^^n f{gtg2) (-^P{a,a+,g,g+)) f^Pia,a+,g,g+)\ (g+gi) 



(190) 



where for simphcity we write P(a,a~^ , g, g^) instead of the full expression 
P{a,a^ ,a* ^a"^* , g, g^) using the notation a, a^ for a, a^,a*, a+*. 

10.2 5.2 Coupled Expansion Coefficients 

The coupled equations for the coefficients specifying the general distribution 
function are obtained from the Fokker-Planck equation (|190l) and as follows. 
For the zeroth order terms 

|PoKa-) = -..(-|-a+^.-)po(«,.-) (191) 
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For the second order terms 

at \ oa oa+ ' 



+ -in (aP^-'^(a,a+)-a+P^''^(a,a+)) (192) 



-P, {a,a ) = -^u.l--a+—a ] P, {a,a ) 



—iujQ P2 ' (a, Oi^) 

{pl'\a,a+)-Po{a,a+)}] (193) 



Hl^ 



-P, ia,a) = -^o.^--a+^-^a )P, (a,a) 



+iujQP2' {a,a ) 

+ -in (a+P2''^{a,a+)-a+P2'^{a,a+)] 

~l'^ (^{^2^'K"+)-^o(a, «+)}") (194) 



-P2 {a,a ) = -^u.^--a+^-^a ) P, (a, a ) 



+—in 10:^^2' (ct, ct^) ^ ctP2 ' (Q^jCt^)) 

:^n(--^P:^•'^{a,a+) + ^P^^■\a,a+) ] (195) 



For the fourth order term 

^ d12;21/ +n • I ^ I ^ + \ d12;21/ +\ 

— P4 ' (a,a^) = -2W ---a + — — a^ P4 ' (a,a^) 

ot \ oa oa'^ } 

-^*"(^^2^'(","+) - ^P2'^'(a,a+) )(196) 

Thus we see that for the general distribution P(Q!,a+, a*, Q!+*,g,g+) the coef- 
ficients satisfy similar sets of coupled equations as in the canonical distribution 
case. 
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10.3 3.3 Rotating Phase Variables and Coefficients 

The transformation to rotating phase variables is as in Eq. (IMl) . The new 
expansion coefficients wih be designated R{/3) where /? — {/3,/3^,/3*,/J~''*} 

and are related to the original -P( a ) via 



Po(^) = Roi^) Pr'^\c^) = Rf^^\A) (197) 

Coupled equations for Rq{P), R^^ ( /? ) and R^ ' ( /? ) are obtained which only 

involve il and A. These are as follows. 
The second order equations are 

d_ 

di 



(^R','\^) - Roi^)) ^+1^^ (^f3R'/{^)~l3+R'/{^)) (198) 



§-^R2''i^) - -^^R'/i^) 



-^in/3 ({i?^^^(^) - Ro{^)} - {i?f (^) - Roi^)}] 
+ izfi(^^{i?^^i(^)-i?o(^)}) (199) 



|i?f(^) = +.Ai?f(^) 



-l^np+ [{rI^H^) ~ Roi^)} {i?r (4) - Ro{^)} 
-^'^[^{R2\^)-Ro{^)}^ (200) 



|(i?f(^)-i?o(4)) = +\^n[(3+R'/{^)~pR',^\^: 

1.^ / d ^1:2/ ^^ , d 0,1 



2^^^-^i?r(4) + ^i?r(4))(2oi) 



The zeroth and fourth order equations are 

d 



g^Roi^)^0 (202) 



^nl2;21/a\ ^nf ^ d1;2/ 



^f (/3) - -:^»« l^^2^'(/3) - T^^f (/3) (203) 
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10.4 3.4 Solution for Distribution Function 

The solution to the coupled equation for the rotating distribution function co- 
efficients is found using the ansatz 

I{ii\^ = exp(/3/?+) $*(/?+) $,(/?) (204) 

where the $*(/?'*') are functions of the /3^ and the '&i(/3) are functions of the /3, 
and where $i(/3) satisfy the coupled equations 

l-^^m = i.A<I>r(/3)4.0(|)a>.(/3) 

^$2(/3) = -izA$2(/3)-izr)/3ci>i(/3) (205) 

This ansatz is consistent with the original equations (I198P - (I20ip for the i?2"' ( /3 ). 
As indicated previously we set i?o(/3)) =0 for the one atom case. 

Following a similar proceedure as for the canonical distribution function we 
find a solution for the $i (/3) that is single valued and not divergent at the origin 
which is of the form 

$i(/3) = /3™ M,„ cos Krat + B,„ sin \v„A (206) 

where 



v„, = \J-m 02 + A2 m = 0,1,2,.. (207) 

However these are not the frequenies associated with population and coherence 
oscillations in the one atom Jaynes-Cummings model. In fact for any detuning 
and Rabi frequency the quantity Vm becomes purely imaginary when ra is large, 
leading to solutions for $i(/3) and ^i(fi) involving hyperbolic functions. For 
example with zero detuning we have 

$i(/?) = /3™ (A„,cosh-^/^nt-iBmSmh-y^nt\ (208) 

These solutions diverge as t becomes large. This indicates that the standard 
distribution function obtained from the Fokker-Planck equation via the standard 
correspondence rules becomes infinite for large time. The corresponding solution 
for ^2{l3) is obtained from Ea. (j205p and thus 

$2(/3) ^^|3^ (^cos -.,^t I 1 + sm -.„,i | 

(209) 
However, a solution with m = leads to a singular /3~ behaviour, so it follows 
that m is restricted to the positive integers. Also, as the ansatz equations are 
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linear the general solution is a sum of terms with differing m so that we finally 
have the solution in the form 



m=l ^ 

A /o^ ■ Sr^ a(m-l) f 1 ^ ( VmBrn + i^Am\ , . f 

*2(/3) = I 2^ Z?^" 'Mcos-z^„tj ^^ '>+sm-iy„it 



n 



Overall, this solution in (J204I) and (I210p for the distribution function is 
unsatisfactory. Not only does it diverge for large t but also the distribution 
function found diverges at all times as /3, (3~^ become large due to the exp(/3/3~'') 
factor. This contradicts the requirement in deriving the Fokker-Planck equation 
that the distribution function goes to zero for large 13, I3~^ . The conclusion is 
that the Fokker-Planck equation (|190l) obtained via the standard correspondence 
rules must be invalid, since if it was valid then the general solution to it that we 
have found would have produced a distribution function that was not divergent 
either in phase space or in time. There is also the question of how to choose 
the quantities An and -B„ so that a valid overall distribution function at time 
i = is obtained. The form of the initial distribution function is not known, 
since the canonical distribution function is not being used in this case where 
the Fokker-Planck equation is based on the standard correspondence rules. Of 
course the only requirement is that the choice of distribution function gives 
the correct characteristic function - the latter being uniquely determined from 
the density operator. Since the characteristic function is obtained via the phase 
space integral of Eq. (|49p it is clear that there is a problem - the bosonic integral 
of the distribution function will not converge due to the exp(/3/3''') factor. No 
matter how the quantities An and -B„ are chosen no valid positive P distribution 
function can therefore be found via the solution given by (J210I) and (|204p . 



(210) 
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